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1 Introduction 



In a previous work p] we proved that the static and axisymmetric solutions of 
the Einstein vacuum equations with a finite number or Relativistic Muhipole 
Moments (RMM) [2] can be described by means of a function u which resembles 
exactly the Newtonian multipole potential if we write the metric in a so-called 
MSA (Multipole Symmetry Adapted) system of coordinates. In fact, this is 
properly the definition of these coordinates, which are constructed iteratively in 
terms of the Weyl coordinates. 

The goodness of these coordinates arises from the fact that they allow us to 
generalize some theorems hold in Newtonian Gravity (NG) that establish the 
existence of certain symmetry groups of the equations satisfied by the classical 
potential whose group-invariant solutions represent the solutions with a finite 
number of RMM in NG The existence of these systems of coordinates is 
proved to be related with the fact that the relativistic equations analogously ad- 
mit those symmetry groups that lead to the Pure Multipole Solutions in General 
Relativity (GR) (those with a finite number of RMM). In addition, the symme- 
try group for the case of spherical symmetry allows us to determine univocally 
the specific MSA system of coordinates by means of a Cauchy problem, whose 
solution provides the standard Schwarzschild radial coordinate. 

These characteristics make the MSA system of coordinates become a very 
relevant reference to describe the gravity of stellar compact bodies whose multi- 
pole structure is known, or suitably estimated a priori, and differs slightly from 
the spherical symmetry. Moreover, a new feature of these coordinates reveals 
more relevance for the description of the Pure Multipole Solutions because they 
provide us with a procedure to establish measurements of high physical and as- 
tronomical interest about the behaviour of test particles orbiting into this kind 
of space-time. 

Until now, only the spherical symmetry solution of Einstein vacuum equa- 
tions has been written in MSA coordinates: this is the Schwarzschild metric in 
standard coordinates. But, do not we have a nonspherical axisymmetric solu- 
tion written in these systems of coordinates, and it would be a very relevant 
success for the description of gravitational effects derived from deviations of the 
spherical symmetry. We are able to write explicitly all the metric components 
of a static solution characterized by any finite number of RMM, in particular 
the monopole, quadrupole and 2^-pole moments, in a system of coordinates 
such that the metric recovers the Schwarzschild limit when all RMM higher 
than monopole vanish. Hence, we are introducing a system of coordinates that 
generalizes the Schwarzschild standard coordinates and it can be used to study, 
in analogy with the spherical case, the physics of test particles in space-time 
slightly different from the Monopole Solution. 

A very useful tool for the study of orbital motions in a classical gravitational 
problem is a second order differential equation called the Binet equation, from 
which the Kepler laws, for instance, can be deduced. From the geodesies of a 
space-time we can obtain a relativistic Binet equation whose resolution allows 
to determine relativistic corrections to Keplerian motion like the correction to 
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perihelion and node line precession. The calculation of the Binet equation in an 
MSA system of coordinates is specially suitable to describe the influence and 
relevance of the different RMM in such corrections. It provides the gravitational 
effects due to deviation from sphericity when comparison with the Schwarzschild 
corrections is made. 

As is known, (see [1] and references therein) the GR theory predicts a cor- 
rection of the Newtonian movement which can be interpreted by means of a 
classical potential type r~^. The Schwarzschild spherically symmetric solution 
of Einstein vacuum equations corresponds with a perturbed Hamiltonian in the 
Kepler problem. In fact, the Binet equation shows directly this relation between 
the orbital equation for geodesies in a relativistic space-time and the orbits ob- 
tained for a classical Hamiltonian dynamical system. 

The question that promptly arises is whether any other solution of the static 
and axisymmetric Einstein vacuum equations can be identified with certain clas- 
sical system by means of a suitable perturbative potential to obtain equivalent 
equations for the orbits. This is one of the purposes of this work which has been 
accomplished through the determination of the equivalent potential associated 
to a perturbed Kepler problem that leads to the same orbital equation as the 
one for certain geodesic of any static and axisymmetric space-time. In this sce- 
nario the study of geodesies of a test particle for a Pure Multipole Solution is 
particularly relevant, because the multipole characteristics of that solution allow 
identification the contributions of its RMM within a classical description of the 
problem. And here is where MSA coordinates become absolutely prescriptive. 

Standard techniques for the study of the equivalent classical problem can be 
performed. In particular, the first order perturbation theory throw interesting 
results about relativistic corrections to the perihelion in Keplerian motion. As 
we will see, expectations arise from this study when the role of the quadrupole, 
and the magnitude of its contribution, are considered, because certain discrep- 
ancies are obtained with respect to other similar calculations [5], [6]. In fact, 
experiments for the measurement of the quadrupole moment can be outlined 
from the contribution of the monopole and quadrupole moments to the perihe- 
lion shift at first order. 

This work is organized as follows. In section 2 a procedure to write any 
Pure Multipole Solution in MSA coordinates is shown and, in particular, all the 
metric components of a static and axisymmetric solution with a finite number 
(Mo, M2, Mi) of RMM are exphcitly obtained. In order to do that, the MSA 
system of coordinates adapted to this kind of solution is previously calculated. 
The good behaviour of the metric and the structure of it, for a general case, is 
discussed. 

In section 3 we calculate the geodesies of such a metric. We show that the 
restriction to the constant equatorial hyper-surface leads to a geodesic equation 
which can be written as a one-dimensional equivalent problem for an effective 
potential. 

Section 4 is devoted to introduce the Binet equation from its classical deriva- 
tion, and we show the relativistic Binet equation for the Pure Multipole Solution. 
We derive the perturbative potential that identifies the calculation of geodesies 
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for any Pure Multipole Solution with the resolution of a classical Hamiltonian 
perturbation of the Kepler-problem. 

In section 5 we show the relevance of the quadrupole moment on the rela- 
tivistic correction to the perihelion for a test particle in a perturbed Keplerian 
orbit. We put our attention on the effect of the quadrupole moment at the 
same order as the monopole correction introduced by the Schwarzschild solu- 
tion. This result suggests the possibility of measurement of the quadrupole 
moment through this contribution. 

Finally, we discuss on the results obtained in a conclusion section. 

2 The Pure Multipole Solutions in the MSA sys- 
tem of coordinates 

The static and axisymmetric solutions of the Einstein vacuum equations with a 
finite number of Relativistic Multipole Moments (RMM) will be referred to as 
the Pure Multipole Solutions from now on. Some authors have devoted works 
on researching about those solutions [2], [8]. 

The MSA (Multipole Symmetry Adapted) systems of coordinates were de- 
fined [T] as those that allow writing the metric component in terms of a function 
u resembling the Newtonian gravitational potential, or equivalently, we can say 
that the expansion of that metric component (ggo = — 1 + 2u) in power series of 
the inverse radial MSA-coordinate provides a multipole Thorne structure with 
vanishing Thome's rests [5] . The MSA system of coordinates belongs to a class 
of ACMC introduced by Thorne [5], with a suitable choice of asymptotical be- 
haviour for the case of equatorial symmetry (the odd order RMM are null) . The 
system of MSA coordinates {i"} = {t,r,y = cos9,(p} {a = 0..3) can be con- 
structed iteratively in terms of the Weyl coordinates {x"} — {t, R,uj = cos 9, ip} 
as follows [1]: 



r^R 



OC ^ 



R" 

n—1 



y = w -\- 
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E^nMi^. (1) 



n=l 



For the purposes of this work we consider the gauge ([T]) up to certain order of 
approximation 0{1/R), and we proceed to calculate the system of coordinates 
associated to a Pure Multipole Solution with the set of RMM desired. In par- 
ticular, the constructive method showed in [1] is followed here to calculate the 
MSA coordinates associated to a solution only possessing the three first Mul- 
tipole Moments of mass: Mq = M, M2 = Q and M4 = D. In the Appendix 
we show explicitly the expressions of the functions fn{^) and gn{^) up to order 
O (l/R^) in the power of the inverse radial Weyl-coordinate R ((58)) - (|59p . The 
solution managed can be considered an exact Pure Multipole solution in the 
sense that it only poses those RMM, but nevertheless the components of the 
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metric are written in terms of a power series of the inverse radial coordinate 
r up to the order of approximation of the gauge ([T]), except for the goo whose 
analytic expression is as follows: 



2 

.900 ^1 + — 



,n=0 



(2) 



where + 1 is the number of RMM of the Pure Multipole Solution selected 
[N = 2 for our case) and P2n (y) stand for the Legendre polynomials depending 
on the MSA angular variable y = cos 6'. 

We first need to perform the inverse change of coordinates ([T]) to write the 
metric in the MSA coordinates from the general Weyl line element of a static 
and axisymmetric vacuum space-time: 

The case of spherical symmetry is especially relevant because we know ex- 
plicitly the sum of the series appearing in the gauge (P) and the standard radial 
coordinate of Schwarzschild is recovered (see [T] for details). The MSA ({r, y})M 
system of coordinates for the spherical symmetry is the following: 

R 

r = M + — (r+ + r_) = Af(a; + 1) 
R 



where r± = vT±2cJA + A^, A = M/R, and {x,yp} are the prolate spheroidal 
coordinates [10] , [11] . The inverse relations between these coordinates are given 
by the following expressions: 



R = rsJl + y^X^ -2\ 
1 + y2A2 - 2A 



(5) 



where A = Al/r. By performing the change of coordinates in the line element ([3]) 
we obviously obtain the known Schwarzschild metric in standard Schwarzschild 
coordinates: 

els'" = - (^1 - dt^ + (^1 - ' dr^ + r^dO^ + sin2 Odip^) . (6) 

The explicit expressions (|3|[S]) are shown because we are able to recover these 
results from the general case of a Pure Multipole Solution by taking all RMM 
equal to zero except for the Monopole, as we will see in what follows. 
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The inversion of the series in the gauge ([T]) to obtain relations of the type 
R = R{r,y), to = uj{r,y), for a general case, is a straightforward but cum- 
bersome calculation. In the Appendix we show the inverse relationfQ ([5(I | -(PT |) 
corresponding to the previously calculated MSA coordinates associated to a 
Pure Multipole solution only possessing the RMM Mq, M2 and M4. 

The transformation of coordinates leads to the following components of the 
metric: 



9ry{x) = 
9yy{^) ^ 

9tt{x) = 
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(7) 



where F(x) = e^2>i'(a:)+27(a:)^ ^YiQ subindices denote the derivative with re- 
spect to that coordinate. 

By using the expression (|59p into ^ we obtain the following metric compo- 
nents, up to order O {^''^ (although we only show here the first terms because 
the expressions are too large) in powers of the inverse MSA radial coordinate r 
(A EE M/r): 



1 2A + + [8 -)- q{l - 3y2)] A^ - y [-14 - 3Mg(l - 3?/^)] A^ 
[-896 + d(2695/ - 2310?/^ + 231) + g(-240 + 72Qy^) + 
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— [d(1715?/^ - 1785y2 + 210) + q{-A2y^ - 94)+ 
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^ [d{89915y^ - 91560^^ + 10605) -I- g(-484?;2 _ 2508)- 



^It can be seen that the expansion of the expressions ((SJl in power series of A reproduces 
this results for the Monopole case. 



6 



+ (-9929 V + 98616?;2 _ 9864)] 
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(10) 



where q = Q/M^, d = D/M^ are dimensionless parameters representing the 
quadrupole and 2'*-pole moments respectively. 
Let us briefly analyze these results: 

i) First, we can see that the cross term of the metric gi2 vanishes up to this 
order of approximation, and in fact, the preservation of the diagonal aspect of 
the metric is a characteristic of the system of MSA coordinates. Supposing that 
the Jacobian of the coordinate transformation is regular and the inversion of 
coordinates is able, i.e., det = \rRyu: — r^y}i\ ^ 0, it can be seen that the metric 
component gry vanishes since we can write the following expression 



Fix) 
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-LBi(r,y)^=j.(£), 



(11) 



where LBi{r,y) = r,^^ V,;()V, () = R^dR{)dR{) + (1 - ^'')d^{)d^{) represents 
the Laplace-Beltrami operator with respect to a 3-dimensional Euclidean metric 
(with axial symmetry) written in Weyl spherical coordinates (see [X, for details), 
and it was seen in [I] that the MSA system of coordinates fulfills the condition 
LBi(r,y) = 0. 

a) Second, we can see from the expressions ((8|)- p0|) that the metric written 
in the MSA system of coordinates leads to the Schwarzschild limit by consider- 
ing equal to zero all RMM higher than the Monopole, as well as that the MSA 
system, given by expressions ([55)) - (ICTl) . recovers the standard Schwarzschild co- 
ordinates. 

Hi) And finally, the 500 [2; = 2;(x)] metric component results to be equal to 
([2|) as claimed by that system of coordinates. 



A general expression for the metric components of the Pure Multipole so- 
lution with a finite number of RMM (Af, Q, D) in MSA coordinates is the 



7 



following: 
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A2 
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(12) 



where Pn{y) stands for the Legendre polynomials, and the Ui, Di and Ti denote 
polynomials in the angular variable y of even order depending on the higher 
RMM {Q and D for this case). 



3 The geodesies of the Pure Multipole solutions 

We proceed now to calculate the geodesies of the metric associated to a Pure 
Multipole solution with mass, quadrupole and 2^-pole moments (M, Q, D re- 
spectively) in the corresponding MSA system of coordinates. The set of equa- 
tions for the geodesies is the following 

dH dtdr^ dtd§^ 

ds^<^''+d-sd-s^^'^'' + d-sd-s^^''' = ' ^ 
d'^r f dr\ ^ ^dr d9 ^ f dt^ 



(14) 




d^e j d9\ ^ ^drde^ ( dtV ^ 

^d?''' + U + ^TsdS^^''' - \js) + 



2 . . 2 



-\$) ^.-^--(S)'^.-^33=0 (15) 

d^w dw dr ^ dw d,Q ^ 

where s denotes the afBne parameter along the geodesic. Since the metric is 

static and axisymmetric, both Killing vectors, ^ and ry, representing the cor- 
responding isometries, remain constant along the geodesies, i.e., ^"Zq, = /i, 
rfza = I, z" being the tangent vector to the geodesic, h, I are constants, and 
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therefore the Eqs. (IT51) and (IT51) become goo-r — h and gzs.-y- = / respectively. 

as as 

The norm of the tangent vector Zj3 = e can be written as follows: 



2 ' 2 



s 

with fc = e , and hence, the geodesic equations ([T4| and (fT5|) can be 

5oo 533 
written as follows: 

a.lnVM+(^^j (^^ja,lngn = ^[9.fc + fca.ln.g22] 

(18) 

fde\\^ , fdr\fde\^^ 1 m , i 

a,lnVM+(^^j (^-j5.1n522 = ^[9,fc + H-lnffn] 

(19) 

Let us consider now the case of geodesies with constant 9 and ^ ^ 0, i.e., 

ds 

they are not radial geodesies but those constrained to a constant hypersurface 
{9 = 9o) with coordinates {t,r,ip}. This restriction is compatible with the Eq. 
(fT^ since it leads to 

[a^fc + /c9^1n5n],-^^„ =0, (20) 

which is fulfilled at the equatorial plane {9n — tt/2 or equivalently y = 0) 
because the expressions (fT2|) allow us to holcijthat ((9y.goo)la=o — {dygii)\y=o = 
{dyg33)\y=Q = 0. Therefore, with respect to the relevant geodesic (IT51) . and by 
using the Eq. (|17p . we obtain on the equatorial plane the following expression: 



ds2 2 51 
This equation can be written as follows: 

|)%V.„ = C, (22) 

Vef f being a so-called effective potential which can be obtained by integration 
as follows 

fc „ , / <7i 1 \ , k 



Veff = I — 9,. In 



V A; / Qii 



511 ^ K / 511 
As can be seen, the Eqs. (P^. (1^^ reproduce the Eq. (IT7)) for constant 



■^Let us note that we have considered equatorial symmetry and hence goo only contains 
even order Legendre polynomials, and the polynomials U, and T depend on even powers of y. 
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4 The Binet equation 



As is known, the classical problem of Kepler consists on determining the move- 
ment of a test particle within a gravitational field generated by a potential of 
the type V{r) n/r. Since this potential is conservative and leads to a gravita- 
tional force orientated towards the center of the source, the conservation of the 
energy E for a particle with mass m and orbital angular moment J is followed 
(/i = GMm): 



1 / • \ 7TI 

E = - { rar'^ + r'^raO'^ + sin^ Oip^ - GM— 

2 V J r 

J = mx A V ^ X ■ J = , (24) 



and the particle moves on a constant plane — 7r/2. Hence, one may define an 
effective potential in the following way 

1 jn 
E=-mr^ + ^,f^ir) , ^ _^ - GAf- , (25) 

and the equations of motion are the following: 

^2^fdrY_2E J2 ^2GM ^ _ / dip\ _ J ^^^^ 

Consequently, the equation of the orbits is given by the following expression in 
terms of a variable u = 1/r: 

du\ 2 2to 



u 



—^{E + GMmu) , (27) 



dip J ,P 

and from it we can easily obtain (by deriving with respect to (p the Eq. (|27|)) a 
second order differential equation for the orbit of the test particle: 

d'^u GMm^ 

This is the Binet equation which can be solved to derive the three laws of 
Kepler concerning the orbit of a test particle describing a closed ellipse around 
the source, for the attractive case E < 0. 

Let us consider now a perturbation of the Newtonian potential of the type 
—a/r^. It is straightforward to see that a generalization of the Binet equation 
is obtained for this case as follows: 

d^u GMrn^ am „ 

Since the corrected potential, i.e., V{r) = — does not depend on the 

angular variables hence, the gravitational force is still central and conservative 
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and we can go on considering the conservation of the orbital angular moment 
of the particle moving along the orbit on a constant surface. 

If we calculate the geodesies of the Schwarzschild metric in standard coor- 
dinates (jH) we obtain, for a constant angular variable 9 = 7r/2, the following 
equation: 



(30) 



where h and I are constants, derived from the isometries as mentioned in pre- 
vious section, that represent the energy and the angular moment per unit mass 
respectively. In fact, this Eq. (150)) is exactly the Eq. (|22p with the effective 
potential for the spherical symmetry case given by the following expression: 

From the Eq. ([50)1 the orbit of the particle is described by the following equation 

— +u = -e- + 3Mu\ (32) 

u = 1/r being the inverse of the standard-Schwarzschild radial coordinate. Since 
I is considered to be the angular moment of the particle per unit mass {I = J/m) 
we can hold by comparing this equation with (j29p the following statement [1]: 
The orbit on a constant surface 9 = t:/2 corresponding to a timelike geodesic 
(e = —I) of the Schwarzschild space-time is given by the same equation as the 
Binet equation corresponding to the Newtonian Kepler problem with a perturbed 
potential of the type —a/r^ with a = j'^M/m. 

As is known [5], [6], the relativistic Binet equation (j32p allow us to calcu- 
late corrections to the Newtonian orbit even for the motion around a spherical 
distribution of mass, which is no longer closed. In particular this relativistic 
effect has been tested in our solar system and it amounts to a slow precession 
of the perihelion of the orbit of Mercury. Predictions of the General Relativity 
correcting the classical gravity, as the precession of the perihelion or deflection 
of light, arises from the resolution of Binet equation and constitute the well 
known ests of GR. 



4.1 The generalized Binet equation 

We want to extend the above-mentioned statement to any static and axisym- 
metric space-time corresponding to a nonspherical mass distribution, by means 
of a generalized relativistic Binet equation. Of course that RMM higher than 
the mass will also contribute to the precession effect and, in principle, these mo- 
ments could be calculated by measuring the precession of a certain number of 
test particles orbiting at conveniently different distances from the gravitational 
source. We will devote last section of this work to discuss this point. 
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Therefore, we need to compare the relativistic Binet equation with a classical 
Binet equation corresponding to certain gravitational potential. According to 
Eq. (12^ we can define an effective potential in such a way that the orbital 
equations (on the equatorial plane) corresponding to geodesies of any static 
an axisymmetric space-times are exactly given by the classical Binet equation 
related to some perturbed Kepler problem. 

Hence, we have established a relationship between the equation for geodesies 
in a relativistic space-time and the orbital equations for an associated Newtonian 
potential. But moreover, since we have constructed Pure Multipole solutions 
with a meaningful physical interpretation in terms of their finite number of 
RMM as successive corrections to the spherical symmetry, the Binet equation 
for these solutions will provide us with the contributions of the different RMM to 
the relativistic effects correcting the Newtonian orbits. In order to support this 
assertion, we must remind that the MSA system of coordinates used to write 
the metric and its geodesies are adapted to the set of RMM of the solution, and 
these coordinates recover the standard Schwarzschild coordinates in the limit 
M2n = 0, Vn > 0. 

In fact, from Eqs. (|22|) and (|23| we have that 



du 
dip 



511 ' 



(33) 



where the notation m = 1/r is used again. Finally, the derivative of the above 
equation with respect to the azimuthal angle ip leads to the following equation: 



d^u 



ld_ 
2du 



533 4 

u 

L511 



g _ ]]^\ 933 ^ 
500/ l'^ 



(34) 



For the case of a Pure Multipole solution with a finite number of RMM (M, Q 
and D) this equation provides the generalized relativistic Binet equation, for the 
orbit on the equatorial plane (y = 0) of a test particle moving along a timelike 
geodesic (e — —1), as follows: 



dip^ 



M 
72" 



3M 



§(5 + 6.') 



30 



15 



D{133 + U0h'^)]u^ + Oiu^) , 



(35) 



where terms in powers of u higher than 5 have been neglected. As can be seen, 
if we take all RMM higher than monopole equal to zero the relativistic Binet 
equation ([32]) is recovered. 

Finally, we want to obtain explicitly which is the perturbed Kepler problem 
associated to this relativistic Binet equation. As we already saw, a perturbative 
potential of the type —a/r^ leads to a classical orbital equation (|29l). that 
reproduces exactly the orbital equation corresponding to a tinielike geodesic 
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of the Schwarzschild metric on the equatorial plane. In analogy with this result, 
if we handle with a Pure Multipole solution with a finite number of RMM, a 
suitable perturbed potential to the Kepler problem can be selected to obtain a 
generalized Binet equation that equals the Eq. (1551) . Nevertheless, despite the 
spherical symmetry case, we will see now that this potential is given up to a 
suitable order of approximation in the power of the variable u. 

Accordingly to Eq. (|25l) let us consider now a generalized effective classical 
potential ^effir) given by this expression: 

*e//-2^-GM^ + V^pM , 1^,(0 = ~^ + y/*"^(r), (36) 

where Vp^*^*'^(r) denotes the new perturbation considered in addition to the 
previously studied. The comparison of the resulting Binet equation with the 
expression (|35|) will supply us with the perturbed potential that provides dif- 
ferent contributions due to the RMM higher than monopole. The conservation 
of energy E and the orbital angular momentum of the particle J lead to the 
following equation for the orbio 

^ J/r^m 

V2/m(£;-$e//) 

or equivalently, we can write the following generalization of the Eq. (j27p : 

+u'^^{E + GMmu + au^^V^^^''"{r)) , (38) 

and therefore, the generalized classical Binet equation is 

(Pu GMm? am o m d o.f»r, , , , 

By comparing this Eq. (15^1) with (I55|) we can hold the following statement: 
Proposition: 

The classical Kepler problem for a mass distribution M , with total energy 
E and orbital angular momentum J , perturbed with a potential given by this 
expression: 

Vp{r) = -au^ ~ QmG{- + 2,h^)u^ + QMmG{ h^)u^ -f 

2m M 8 2 



+ (^-y/iVM^GQ 



(40) 



^Let us note that I I37I1 can be obtained from Eqs. I I26II 
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provides ( up to the desired order in the power of the variable u ) the same orbital 
equation for a test particle of mass m, as the corresponding orbital equation for 
the timelike geodesic on the equatorial plane of the relativistic Pure Multipole 
solution with a finite number of RMM if the following values of the parameters 
are considered: 



1 = 



J 



2E 



1 



(41) 



Let us note that the above expressed potential can be written in terms of 
the dimensionless parameter A = Mu = M/r, which controls the perturbational 
character of it, as follows: 



14 



q Uh^ -6 + 7 



J' 



h Aa\(-3 + 25/1^)+ 
14 ^ ' 

28q^{l + h^)--d{19 + 20h 



(42) 



where q = Q/M^ and d = D/M^ denote for dimensionless parameters associated 
to the quadrupole Q and 2'*— pole moment D respectively. 

Consequently, the Binet equation is given by the following expression: 



+ . = M(!^)^G+[(f +9/.^)(^)^GQ + 

2 



,,^^(mYI 18 150. 



^-^Q+['^] G(30(l + /i2)^+M2g^(-3 



J 



M 



.15. 

y 



^^D{2Qh^ 



19) 



3M 



22h^ 



(43) 



5 Measuring the quadrupole moment 

The Eq. PO)) shows that a quadrupole contribution arises at order , the same 
order as the known Schwarzschild correction due to the monopole. The role of 
the quadrupole becomes relevant enough when the relativistic correction to the 
perihelion of a test particle orbit is calculated. Hence, let us consider the above 
potential (l40l) and let us develop up to first order of perturbation theory the 
following perturbed kepler problem (for 71 = 3): 

V ^-GMm}--^ , H = H^+iHi , Hi = -^, (44) 

where ^ is a small parameter, and Hi is the perturbation of the Hamiltonian 
H. 
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According to the standard theory of perturbations ^ within angular-action 
variables, the averaged rate of secular precession due to the perturbation is given 
by the following expression: 



- 1 

<; = - 

T 



(45) 



T being the average time interval considered, in particular the period of the 
nonperturbed orbit; stands for the angular position of the periastron on the 
orbit plane, and Hi is the time-averaged perturbed Hamiltonian which can be 
calculated by using the conservation of the angular moment (J = mr^ip—cie) 
as follows: 



1 1 m 

—dt = 

r" T J 



7 dip, 



(46) 



where r can be expressed in terms of if by means of the nonperturbed orbit 

1 1 



P 



-(1 -I- ecosip), 



(47) 



e being the eccentricity and p — 



is the so-called ellipse's parameter. 



Therefore, the averaged rate of precession is given (for the n = 3 case) by 



T J* T 



J2 



J4 



(48) 



The first term of (|48l) is the monopole contribution already known, which is 
provided by the Schwarzschild solution, whereas the second term is the correc- 
tion to the precession due to the quadrupole moment. If we take into account the 
value of h (|4T|) and consider J and E in terms of the eccentricity e and the semi- 
major axis a of the orbital ellipse, i.e., = a{l — e'^)GMm? , \E\ — GMm/{2a), 
we have that 

1 Ar 



6tt 

T 



(49) 



where 



M 



a(l 



is a dimensionless parametei|f| less than 1, the quadrupole 



contribution is of order ^■^ and its sign depends on the relative value between 
M and a: 



1 



M 



< 
> 



M « a 
6M > a. 



(50) 



In [5] , |6] the authors look for patterns of regularity in the sign of the contri- 
butions for different RMM, ending up the conclusion that the linear quadrupole 
term is negative for a positive quadrupole. This pattern is verified by our cal- 
culation except for the case that 6M > a, i.e., a test particle closely orbiting 



*Let us note that M must be considered in length units (M ~^ )- case of the 

Sun the value of this length is about 1.476 km, and hence M << a for any planet of the solar 
system. 
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around a strong gravitational source. The discrepancy arises because in [B] the 

authors do not obtaiiH, up to this order in the parameter C, the relativistic 
3M 

contribution appearing in (|49l) and the contribution of the quadrupole at 

a 

order C,'^ is given only by the Newtonian term —1/2. 

Since the first contribution of the quadrupole moment to the perihelion shift 
appears at the same order in A as the monopole contribution, this calculation 
provides a procedure for measurement of the quadrupole moment. In principle 
a test particle conveniently far away (A < 1) from the source should be affected 
only by the monopole and quadrupole contributions which are the only effective 
corrections at that distance (contributions higher than A'^, in the first order of 
perturbation theory, are considered negligible). Dropping any other external 
effects, for a suitable isolated source-particle system acting on the perihelion 
precession of the orbit, the quadrupole contribution can be estimated to obtain a 
tentative measurement of the quadrupole moment since the monopole correction 
is well known. In fact, the relative value between both corrections at this order 
of perturbation theory is the following: 

A = -fri — = ^ ° — ^ = Cg 3 . 51) 

Let us calculate an estimate of this relative value for the orbit of Mercury. 
As is known, Einstein's theory of gravity leads to a relative correction of the 
Newtonian value for Mercury's secular perihelion drift of 42.98"/century [12]. 
The rate of precession for the relativistic monopole (Schwarzschild correction) 
can be calculated from the first term of our Eq. gSl) to obtairl!|38.01"/century 

The quadrupole contribution can be calculated from the second term of 
the Eq. (|49p . and hence, the key point of its estimate, as well as that of the 
relative value (jSip . depends on the dimensionless parameter q related with the 
quadrupole moment of the Sun. This is a matter of increasing researclQ, because 
the quadrupole moment is a relevant quantity for a lot of related measurements; 
in fact, the evaluation of the solar quadrupole moment still faces some contro- 
versy: on one side, the theoretical values strongly depend on the solar model 
used, whereas accurate measurements are very difficult to obtain from observa- 
tions, in particular the value of the quadrupole moment can be inferred to be in 
agreement with the experiment observations of precisely the perihelion advance 
of Mercury (see [M] and references therein). 



^Let us note that our parameter ^ is exactly equal to the parameter used in [6] to develop 
the expansion series. 

^The quantity {6tt/t)C, provides an estimate of the Schwarzschild correction, with the 
following values of the parameters: eccentricity of the orbit e = 0.2056, the semi-major axis 
of the orbit a = 6.4529 X lO^fcm, M being a half of the radius of Schwarzschild for the Sun 
M = GMq/c'^, and period of the Mercury's orbit t = 7600428s. 

^The gravitational multipole moments describe deviations from a purely spherical mass 
distribution. Thus, their precise determination gives indications on the solar internal structure. 
It is difficult to compute these parameters and the way to derive the best values and how they 
will be determined in a near future by means of space experiments is the aim of several works 
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Different estimates of the quadrupole moment of the Sun [T3] , [13] , provide a 
range of values from a theoretical determination J2 = (2±0.4) x 10"'' to another 
bounds around J2 = (1.4±1.5) x 10~^ [H]; nevertheless, the quadrupole moment 
of the Sun may not exceed the critical value of J2 = 3.0 x 10~^ according to 
the argument given in [15], based upon the accurate knowledge of the Moon's 
physical librations (these data reach accuracies at the milliarcsecond level). 

The factor M/a for the Mercury's orbit is about 2.257 x 10~*, and the 
parameter C, strongly depends on that value, since 1/(1 — e^) is around 1.04415, 
and so C — 2.356 x 10^®. Therefore, Eq. ([5T|) for the relative value of the 
quadrupole contribution with respect to the Schwarzschild correction leads to 
the following estimate: 

|A| ~ 1.18 X lO^^g , (52) 

or equivalently, the rate of precession due to the quadrupole contribution is 
estimated as follows: 



Gtt 



^'^'-2 + ^T 



4.48 X IQ-^" /century. (53) 



This result is perfectly in agreement with the current predictions of measure- 
ments and values of the advance of Mercury's perihelion deduced from observa- 
tional data: As is known [14], once correcting for the perturbation due to the 
general precession of the equinoxes and for the perturbation due to other planets, 
the advance of the perihelion of Mercury is a combination of a purely relativistic 
effect and a contribution from the Sun's quadrupole moment. One may compute 
the corrective factor to the prediction due to General Relativity (42.98"/century, 
the Schwarzschild contribution, which does not include the quadrupole correc- 
tion): this factor, i.e., the relative value of the solar quadrupole correction with 
respect to the Schwarzschild contribution is 

2.821 X 10"'^(2.0± 0.4) , (54) 

for a value of the quadrupole moment J2 — 2.0 x 10^^, or equivalently, the 
correction to the perihelion precession due to the quadrupole moment is about 

2.425 X IQ-^" / century , (55) 

These calculations are done within the Parameterized Post-Newtonian (P.P.N.) 
formalism describing a fully conservative relativistic theory of gravitation (see 
Eq. (1) in fT4]), and they are in accordance with our Eqs. ([52]) and ([53]) respec- 
tively by taking into account the appropriated relation between the parameter 
J2 in P.P.N, formalism and the relativistic multipole moment Q defined by Ge- 
roch [2], i.e., \Q\ = Mi?| J2. Hence, from ([52]) we have that 

|A| ~ 4.45 X IOV2 ~ 8.9 X lO"'^ , (56) 

where a value of J2 = 2.0 x 10^^ has been used, or equivalently, Eq. ([53]) leads to 
the predicted relations between observational or theoretical quadrupole moment 
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and the quadrupole correction to the perihelion precession of Mercury's orbit: 



iQ ~ 4.48 X lO^Vz - 1-12 X 10^ J2"/ century ~ 2.24 x IQ-^" /century. 

(57) 

6 Conclusion 



In this work we have tried to show the relevance of the MSA system of coor- 
dinates for the description of static and axisymmetric vacuum solutions with a 
finite number of RMM, particularly for those which are slightly different with 
respect to the spherically symmetric solution. In [16] the authors study the be- 
haviour of different geodesies of quasispherical spacetime, for example the case 
of self-gravitating sources with the exterior gravitational field of the M — Q^^^ 
solution [5]. 

First, we have explicitly written a static and axisymmetric vacuum solution 
with a finite number (M, Q, D) of RMM in MSA coordinates. The expressions 
obtained for the metric components are approximated because the MSA system 
of coordinates are constructed iteratively by means of a power series. But, these 
results allow us to handle with the Pure Multipole solutions as generalizations of 
the Schwarzschild solution in the sense that each one of the metric components 
are written as a series whose first term represents the monopole solution and 
the following terms provide the successive corrections to the spherical symmetry 
due to the other multipoles (|12|) . Furthermore, the goo metric component is 
calculated to any order and it resembles, as it was desired, the formal expression 
of the classical multipole potential. Since these solutions are static, we are 
actually giving the Ernst potential |17) of the solution in such a way. Until 
now, no other solution has been written in MSA coordinates except for the 
Schwarzschild solution. The expressions obtained for the metric are supported 
by the calculation of the MSA coordinates in terms of the Weyl coordinates 
(as well as the inverse relations). We are providing a system of coordinates 
that generalizes the standard Schwarzschild coordinates for the cases of Pure 
Multipole solutions. 

Second, we have used this system of coordinates to study the behaviour of 
test particles orbiting in a space-time described by a Pure Multipole solution. 
Two results seems to be specially relevant: One of them is the possibility of 
finding an equivalent classical problem; we are able to write the orbital equation 
associated to geodesies in the equatorial plane identically equivalent to the Binet 
equation obtained from a classical perturbed Kepler problem. In other words, we 
can describe the orbital motion of a test particle in the same way as it is studied 
the classical field equations for a problem of a conservative potential endowed 
with certain perturbation. We calculate explicitly this perturbative potential in 
terms of the physical parameters of the virtual Hamiltonian dynamical system 
(energy E, particle mass m, orbital angular momentum J...) and the RMM of 
the vacuum solution. 
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This relationship between the resolution of a classical dynamical system and 
a geometrical description of the problem by means of an associated riemannian 
metric is an analogous result to the prescription derived from the Maupertuis- 
Jacobi principle |18) : the trajectories of a mechanical system, with a natural 
Lagrangian function L = [l/2)gijq^q^ — V{q), are geodesies of the Jacobi metric 

= 2{E — V)gij for a fixed total energy E of the system. Possibly be some 
connection between those schemes can be explored. 

The other important result of this study of the geodesies in MSA coordinates 
is the calculation of high order relativistic corrections to Keplerian motion. In 
particular, we have calculated the perihelion shift due to the correction of the 
quadrupole moment in addition to the Schwarzschild contribution, and it is 
shown that both corrections arises at the same order in the perturbed potential. 
Classical techniques of perturbation theory can be used to make this calcula- 
tion since the description of the relativistic motion is made from the equivalent 
Newtonian problem. In [S], [B] the authors develop a calculation of relativistic 
corrections to Keplerian motion; the advantage of our work is that we do not 
need to introduce a parameter ad hoc to perform the expansion series, because 
the MSA system of coordinates itself leads to geodesies of the solution written 
in such a way that the corrections due to any RMM are clearly distinguished. 

The estimate of the quadrupole contribution leads to the conclusion that it 
is small compared with the Schwarzschild correction for the case of Mercury's 
orbit, but it is perfectly detectable by present experiments, since one of the 
techniques used for measuring the quadrupole moment of the Sun is just by 
means of the perihelion precession of the orbit [13] , [M] . 

Comparison between both corrections ([5T|) leads to the conclusion that the 
quadrupole contribution to the perihelion precession may be quantitatively sig- 
nificant with respect to the Schwarzschild correction for scenarios with a very 
massive source and test particles closely orbiting around: Eq. p9)) shows that 
it is important not only the value of the quadrupole moment but also the factor 
M/a. So, if we think about an experiment or astrophysical system with a test 
particle closely orbiting around a strong gravitational source we could handle 
with a value of the factor M/a which may compensate the factor q whose value 
is rather small for standard astrophysical systems (except for compact binary 
systems). 

A future generalization of our results to a stationary and axisymmetric case 
will provide us with a more realistic scenarios, and the estimates will become 
very relevant if the calculation is applied for instance to an axisymmetric spin- 
ning star. For example, neutron stars are extremely compact objects typically 
around IMq or 2Mq compressed in a radius of a few kilometers, and hence, the 
corresponding value of M/a for a test particle orbiting closely around a neutron 
star can be near to 10^^. 

In addition, some works have studied the central role of the innermost stable 
circular orbit (ISCO) in the relativistic precession of orbits around neutron stars 
[To] . |20| . Strange stars has also been considered as relevant sources (of QPOs for 
example [21]); Strange stars are objects with two main characteristics: they are 
made of a mixture of quarks and they have no minimum mass. It has also been 
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shown in [Tn],[2I] that ISCO is located outside the strange star for relatively 
high mass (ranged from 1.4 to 1.6 solar masses depending on the different EoS) 
even at very high rotation rates, this fact being the main difference with respect 
to neutron stars. In ^2] Shibata and Sasaki calculated analytically the ISCO of 
neutral test particles around a massive, rotating and deformed source in vacuum, 
including the first four gravitational multipole moments (following the scheme 
given by Shibata the authors in [23] extend it to the electrovacuum case.); it 
is possible to study the role of the quadrupole moment of mass related to the 
neutron stars, whose value is not constrained to be small, and the factor M/a 
can be considered for hypothetical orbits nearby the ISCO (r « 6M). 



7 Appendix 

The gauge ([1]) is defined in terms of series, up to order 0(l/i?^), with the 
following functions /„(aj) and g„(w): 
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16 16 16 16 
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f ( \ - 91 1^ 52341 2 138851 13553439 4 892359 A 
~ M [~520~'^ ^ 24024 40040 ^ 3080 j 

„^ / 84665 6 646781 19850195 4 280173 2 675 c\ 

+ DO [ H 1 u)^ + 

^ V 176 64064 32032 1456 64 J 

^,,5/ 383 4897 , 4825 . 1595 . 715 o\ 

QM^ uj* H oj^ H uj'^ + 

^ V 1056 112 462 28 32 y 

,,0 / 5 315 4 35 2 231 . 429 

cj^ H H 

V 128 64 32 32 128 



+ 



' 128 64 32 32 128 J 
,,-1^/^ 445 3890 2 3171 . 159915 4 273 

+ -C + — TTT^ + TTir^ TT^ 



M^£) \ Lo^ H L 

V2288 143 22 1144 16 

,^2^2/^ 936731 1143841 2 34285679 4 4121 (, 715 « 
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(58) 
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= (1 
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2. 1 

^ lOM 

' ^3696 



[DM(-105a;2 + 45) + QM^(-30a;^ + 8) + Q'^{96lj'^ - 36)] 

) [M^(7623w4 - 6930cj2 + 1155)+ 

+ QM^ (254100;^^ + 7546 - 34188^^) + L»M(9702w'' + 27510 - 69580^^)+ 
+ Q2 (6930a;4 + 52356w^ - 20694)] 

£,7(0;) = (1 - cj2)cj(-_i_) [L»QM(-64680a;'^ + 54600a;2 - 8400) + 
490Af^ 

+ Q^M^{-198A5uj^ + 7024^2 + 957) + QM'^(4410w4 - 3038^^ + 336)+ 
+ M'^L>(25725w'' - 13475^^ - 210) + Q^(34776cj'' - 274680;^ + 4032)] 

gs{oj) = (1 - '^^)w— 4:— [Q3(417624012lj4 - 327877128^^ ^ 48074796)+ 

+ i:>gM(-20270250cj6 - 6934063500;* + 595505850cj2 - 92457090)+ 

+ QM^ (-58558500W® + 100720620w'^ + 4832100 - 45945900^^)+ 
+ (-107357250;^ + 15030015a;* - 5780775w^ + 525525)+ 

+ M*Z)(-37837800a;*' + 226746520a;* - 1371678000;^ + 12251400)+ 

+ M^Q^ (-49549500a;'' - 83972460a;* + 63989772a;2 - 4272228)] (59) 

The following expressions show the relation of Weyl coordinates {.R, a;} in 
terms of the MSA coordinates {r,y} up to order 0(A®): 
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